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Use of Quadrupolar Nuclei for Quantum Information processing by Nuclear Magnetic 
Resonance: Implementation of a Quantum Algorithm. 

Ranabir Das^ and Anil Kumar t>+ 

t Department of Physics, * Sophisticated Instruments Facility 
Indian Institute of Science, Bangalore 560012 India 

Physical implementation of Quantum Information Processing (QIP) by liquid-state Nuclear Mag- 
netic Resonance (NMR), using weakly coupled spin-1/2 nuclei of a molecule, is well established. 
Nuclei with spin>l/2 oriented in liquid crystalline matrices is another possibility. Such systems 
have multiple qubits per nuclei and large quadrupolar couplings resulting in well separated lines 
in the spectrum. So far, creation of pseudopure states and logic gates have been demonstrated in 
such systems using transition selective radio-frequency pulses. In this paper we report two novel 
developments. First, we implement a quantum algorithm which needs coherent superposition of 
states. Second, we use evolution under quadrupolar coupling to implement multi qubit gates. We 
implement Deutsch-Jozsa algorithm on a spin-3/2 (2 qubit) system. The controlled-not operation 
needed to implement this algorithm has been implemented here by evolution under the quadrupolar 
Hamiltonian. This method has been implemented for the first time in quadrupolar systems. Since 
the quadrupolar coupling is several orders of magnitude greater than the coupling in weakly coupled 
spin-1/2 nuclei, the gate time decreases, increasing the clock speed of the quantum computer. 

I. INTRODUCTION 



^ ■ In 1985 David Deutsch suggested an algorithm which demonstrated the use of "massive quantum 
! parallelism" inherent in quantum systems [1]. Better known as Deutsch-Jozsa (DJ) algorithm, it 
can distinguish between a 'constant' and a 'balanced' function in an N-qubit system in one function 
call, where as its classical counterpart requires on the average (2 N ~ l + 1) function calls [1,2]. Over 
the years quantum information processing have been demonstrated in various physical systems [3,4]. 
Nuclear magnetic resonance (NMR) has also successfully demonstrated several avenues of quantum 
information processing [7-25]. Quantum algorithms like Deutsch-Jozsa algorithm, Grover's search 
algorithm and Shor's prime factorization algorithm have been successfully implemented by liquid 
state NMR using molecules having weakly coupled spin-1/2 nuclei [10-18]. In such systems each 
nucleus is identified as a qubit and the coupling between the qubits (nuclei) is mediated through 
covalent bonds (indirect spin-spin J-coupling). 

A growing appreciation among researcher's is the use of quadrupolar nuclei with spin>l/2 as a 
suitable candidate for quantum information processing [20-25]. The energy levels of a quadrupolar 
nucleus are equispaced in a liquid, yielding degenerate single quantum NMR transitions. This 



degeneracy is lifted in a liquid crystalline matrix yielding 21 well resolved transitions, allowing the 
21 + 1 eigenstates of a half integer spin I nucleus to be treated as states of an N-qubit system, 
provided (21 + 1) = 2^. In such cases a single quadrupolar nucleus acts as several qubits [21]. 
In such systems, while the quadrupolar splittings are of the order of several KHz, the line widths 
are only of few Hertz. Short and precise transition selective pulses can be applied to such systems 
[21,24]. 

The Hamiltonian of a quadrupolar nucleus partially oriented in liquid crystalline matrix, in the 
presence of a large magnetic field B and having a first-order quadrupolar coupling is given by [26] 

n = n z + n Q = -u i z + J^ Q _ i} (3i 2 - P)s 

= -c / 2 + A(3/ 2 -/ 2 ), (1) 

where ujq = ^Bq is the resonance frequency, 7 being the gyromagnetic ratio, S is the order param- 
eter at the site of the nucleus, e 2 qQ is the quadrupolar coupling and A = e 2 qQS/(AI(2I — 1)) is 
the effective quadrupolar coupling. Though e 2 qQ is of the order of several MHz, a small value for 
the order parameter (S) converts the effective quadrupolar coupling 'A' into several kHz. Prepara- 
tion of pseudopure states, implementation of logic gates and half-adder /subtracter operations and 
quantum simulations have already been demonstrated in such systems [21-25]. However, so far 
only logical operations which do not require coherent superposition have been implemented in such 
systems. In this work, we demonstrate that such systems can also be utilized for quantum infor- 
mation processing by implementing algorithms which need coherent superpositions of states such as 
Deutsch-Jozsa algorithm. Moreover, we propose the use of evolution under quadrupolar interaction 
for implementation of such algorithms. The Hamiltonian of Eq.(l) has two parts- (i) Zeeman part 
(ujqI z ) and (ii) scaled quadrupolar part (A(3/ 2 — J 2 )). The pulse sequence t/2 — (tt)—t/2 focuses the 
Zeeman interaction but allows the the system to evolve under the quadrupolar interaction. Similar 
to J-coupling, quadrupolar coupling provides interaction among multiple qubits, and can be used by 
such sequences to implement multi qubit gates. Since the gate time is inversely proportional to the 
strength of the interaction, and the scaled quadrupolar coupling is three order of magnitudes greater 
than J-coupling, the gate time decreases thereby increasing the clock speed of the quantum computer. 
However, in such systems the relaxation times are also smaller by two orders. Thus decoherence 
takes away some of the advantage of faster gate speeds. We have experimentally implemented DJ 
algorithm using quadrupolar coupling in 23 Na (spin-3/2) nuclei. 



II. DEUTSCH-JOZSA ALGORITHM 



The DJ algorithm determines the type of an unknown function when it is either constant or 
balanced. In the simplest case, f(x) maps a single bit to a single bit. The function is called constant 
if f(x) is independent of x and it is balanced if f(x) is zero for one value of x and unity for the 
other value. For N qubit system, f(xi,x 2 , ...xjv) is constant if it is independent of X; L and balanced 
if it is zero for half the values of Xi and unity for the other half. Classically it requires (2 7V_1 + 1) 
function calls to check if f(x\,x 2 , ...xn) is constant or balanced. However the DJ algorithm would 
require only a single function call [1,2]. The Cleve version of DJ algorithm implemented by using a 
unitary transformation by the propagator Uf while adding an extra qubit, is given by [27], 

\x 1 ,x 2 , ...x N ) \x N+1 ) — ► \x u x 2 , ...x N ) \x N+1 © f{X!, x 2 , ...x N )) (2) 

The four possible functions for the single-bit DJ algorithm are given in Table 1. 



TABLE 1. Constant and Balanced one-qubit functions. 





Constant 


Balanced 


X 


h h 


fs fi 





1 


1 


1 


1 


1 



The unitary transformations corresponding to the four possible propagators U f are 
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(3) 



In the case of 2 qubits there are 2 constant and 6 balanced functions. For higher qubits the functions 
are easy to evaluate using Eq. [2] . 



III. DJ ALGORITHM IN SPIN-3/2 SYSTEM 



The Cleve version of the algorithm implemented here requires two qubits, an input qubit and a work 
qubit [27]. Previous workers have demonstrated the algorithm using two weakly coupled spin-1/2 
nuclei [10,11,13]. Here we implement it on a spin-3/2 nucleus. The energy level diagram of a spin-3/2 
nucleus corresponding to the Hamiltonian of Eq. [4] is shown in Fig. 1(a). The four energy levels 
are labeled as levels of a two qubit system (Fig. 1(a)). There are three single quantum transitions 
among the four levels, labeled as c^oo-oi) ^oi-n an d cj 11 _ 10 , Fig. 1(b) (the subscript denotes the 
energy levels between whom the transition takes place). While the above three transitions are single 
quantum transitions (Am = ±1) they are also single qubit flip transitions. However, according to the 
present labeling scheme the single qubit flip corresponding to 1 00) <-> 1 10) is a forbidden Am = ±3 
transition. The above labeling scheme was chosen to optimize the experimental implementation of 
the algorithm. As demonstrated in 1=7/2 systems elsewhere, optimum labelling schemes can be 
chosen for given logical operations [25]. 

The quantum circuit for implementation of DJ algorithm followed in this work (Fig. 2), is similar 
to those used by previous workers [10,11]. The first qubit is the input qubit whereas second qubit is 
the work qubit. The algorithm starts from a pure state = 1 00) followed by a Hadamard transform 
[11]. A pseudo-Hadamard operation can be implemented by using a high power, low duration 'hard' 
(n/2) pulse along (-y) axis which creates superposition of all the qubits [10,11]. In the spin-3/2 
system the operator of hard {it /2)- y pulse is of the form; 

/ 1 



exp(H y n/2) 



1 \ 

i / 



(4) 



1 -y/3 -1 1 

2V2 v/3 -1 -1 

V -1 V3 -y/3 
where the y-component of spin angular momentum in spin-3/2 system is 

/ -y/3/2 \ 

VS/2 0-10 
1 — s/3/2 

V y/3/2 J 

The state of the system after (n/2)- y pulse is = exp{il y n /2)\ip) = ^[|00) - y/3\01) + y/3\ll) - 

|10>] = 5^[(|0>-|1»(|0>-V3|1»]. 

It is to be noted that unlike weakly coupled spin-1/2 nuclei, the operator of (ir/2)_ y pulse used here 
does not create uniform superposition. However, it does create a coherent superposition of all the 



(5) 



states which can be utilized for 'quantum parallelism' as desired by the algorithm. This also works 
for higher spin systems like spin-7/2 nuclei which can act as 3-qubit system [28]. After creation of 
\ip') we apply the unitary operator Uj which yields \ip'j)—Uf\ip'). The operator Lq is unity operator, 
yielding \ifj'{)=U l \ifj') = ^[(|0) - |1))(|0) - v^l))]. The operator U 2 flips the state of the second 
qubit, yielding; \^) = U 2 \tp') = 573 [(|0) - |l))(-v / 3|0) + |1))]. U 3 flips the state of second qubit 
only when the state of first qubit is |1) and U4 flips the state of second qubit only when the state 
of first qubit is |0) (Table 2). The operators U3 and U4 are thus controlled-NOT gates. It may be 
noted that similar to spin-1/2 case, the information about the function is encoded in the relative 
phase of the two states of the input qubit; (-1) for constant and (+1) for balanced functions (Table 
2). 

TABLE 2. The function (/), the operator Uf, and the wave function \^"f) for one qubit DJ. 
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Density matrices of the system confirms the different functions. The density matrices correspond- 
ing to the states M) are a'L given by 
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The signs of input qubit coherences |11) <-> 1 01) and 1 10) <-> 1 00) , are negative for a" and c^' but 
positive for a' 3 ' and 04, indicating respectively constant and balanced functions. 



After Uf, one needs to make a measurement (Fig. 2). Theoretically this step needs a Hadamard 
gate followed by a readout of input qubit. In NMR, the Hadamard is replaced by a pseudo Hadamard, 
which can be implemented by a (tt/2) pulse. Similarly, the readout is also another (n/2) pulse. These 
two pulses cancel each other and hence in NMR, the result of D J algorithm is directly available after 
implementation of Uf [11]. As seen from Eq. [6], in the signal acquired immediately following Uf] 
the resonance of input qubit at cj i-ii (the central transition of Fig.l) will be of the same sign 
as the resonances of work qubit at cj o-oi an d f^n-io (the outer transitions of Fig.l) for constant 
functions, and of opposite sign for the balanced functions. It may be mentioned that only one of the 
transitions of the input qubit namely 1 01) <-> |11) is observed here, the other transition 1 00) <-> 1 10) 
being Am = ±3. 



The DJ algorithm is experimentally implemented here on 23 Na (spin-3/2) nuclei of a lyotropic 
liquid crystal composed of 37.9% sodium dodecyl sulfate, 6.7% decanol, and 55.4% water. The 
liquid crystal had a nematic phase at 299K [21,31]. All experiments were performed on a DRX 
500MHz spectrometer. Fig. 1(b) shows the equilibrium spectrum consisting of three lines, with a 
effective quadrupolar coupling (A) of about 16 kHz and integrated intensity ratio 3:4:3. 

The 1 00) pseudopure state is created by applying a selective population inversion (ir) pulse on the 
1 10) <-> 1 1 1 ) transition followed by a population equilibration (tt/2) pulse on 1 01 ) <-> |11) transition 
and a gradient pulse to kill created coherences [9,24]. Transition selective pulses are long duration, 
low power r.f pulses applied at the resonant frequency between two energy levels, which excite a 
selected transition of the spectrum and leave the others unperturbed. Let us consider a two level 
sub-system \i) and \j), whose equilibrium deviation populations are and pj 



The operator of a transition selective (9) pulse about y-axis between these two levels would be 



IV. EXPERIMENT 
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A population equilibration (vr/2) pulse will equilibrate the populations and create coherences of the 
form; 

1/1 \\ (Pi §\ \ (\ -V 
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= /(Pi+^-)/2 fe-p0/2\ 
\{pj-Pi)/2 (pi+Pi)/ 1 !)' 
which followed by gradient will retain the populations but destroy the coherences. 
After the creation of pseudopure state the coherent superposition of both the qubits are created 
by a non-selective (7r/2)_ ?/ pulse. At this stage one can apply the various Uf. The function U\ needs 
no pulse and the result given in Fig. 3 indicates that all the three transitions are of same sign and 
hence it is a constant function. The operator f/ 2 in Eq. [3] requires two transition selective pulses 
(^/v^)r >H0 V/v^)l 10>Hll> , where, 
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Here we have used the fact that 
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The result given in Fig. 3 confirms that / 2 is also a constant function. While implementing U 3 of 
Eq. [3] we note that the Pound- Overhauser CNOT gate [8] is similar to f/ 3 , but differs from its exact 
form by a controlled phase operator 
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The Pound-Overhauser CNOT gate is implemented by a transition selective (7r/v / 3)'^^' 11 ^ pulse. 
The controlled phase shift operator of Eq. [13] can be realized by using (i) only transition selective 
pulses or (ii) transition selective pulses along with a evolution under quadrupolar coupling. 

(i) Transition selective pulse method: Transition selective z-pulses can be used to introduce 
specific phases to the different states [30]. For example, the introduces a phase shift of 20 be- 

tween the states 1 01) and |11). is implemented using three selective pulses (it / A) y (4>) x ('n / A)_ y 

on the transition |01) <-> |11); 
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where the x and y-component of the operator of transition between the states |01) <-> |11) are 
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Hence the controlled phase shift operator of Eq.[13] can be achieved by using a cascade of three 
transition selective z-pulses 
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(ii) Evolution under quadrupolar coupling method: Similar to J-coupling, evolution under 
the quadrupolar coupling rotates the system about the z-axis, introducing specific phases to the 
states. The quadrupolar Hamiltonian in spin-3/2 system is of the form 
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(17) 



The operator corresponding to evolution under the quadrupolar Hamiltonian for a time r is 
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(18) 



Hence the controlled phase shift operator of Eq.[13] can be achieved by a combination 
(e lHQT )(n/2) z 1 ^ 11 , where r = 7r/12A is the time period of evolution under the quadrupolar Hamil- 
tonian. 

We have implemented approach (ii) in our experiments. This is because, the use of evolution under 
quadrupolar coupling reduces the number of transition selective pulses, enabling fast computation 
and less errors due to relaxation. All the experiments were carried out with the carrier frequency of 
the r.f. pulses matching with the central transition (on- resonance) . In this situation, the evolution 
in the rotating frame takes place only under quadrupolar Hamiltonian and the Zeeman term does 
not evolve. U3 was implemented by a pulse sequence 
(7r/v / 3)L 1 ^ |11> (e^r)( 7r/2 )|oi>«|ii> 
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(19) 



where r = 7r/12A. The operator U4 was implemented by a similar pulse sequence 
(vr/v / 3)L°° H|01> (e m a r )(7r/2)L°° H|01> , with the same value of r. The result after applying U 3 and 
C/4 are given in Fig. 3, in which it is seen that the sign of central transition is opposite to that of 
the outer transitions, indicating that and f± are balanced functions (Eq. [6]). 

The selective excitation in this paper is achieved using Gaussian soft pulses [32] of length 123/xs. 
During the selective pulses the unexcited transitions continue to experience quadrupolar interaction, 
resulting in a rotation around the z-axis, which leads to phase errors. To minimize such errors, 
the length of the selective pulses were so chosen that the phase rotation is in multiples of 2n [11]. 
However, errors due to relaxation could not be avoided. For example, the peak intensities are 
slightly different from the expected. It is because the relaxation times are: Ti = 16ms for all the 
three transitions; T2 was 14ms for the central transition and 4ms for the outer transitions. Since, the 
relaxation time of the outer transitions is less than the central, the coherences of the outer transitions 
decay faster, decreasing their peak intensities. 

In quantum information processing by NMR, the gate time is of the order of the inverse of coupling 
and the coherence time is proportional to inverse of the linewidth. In liquids the coupling values 
are ~ 100 Hz and the linewidth ~ 1 Hz yielding a gate time ~ 3 ms, a coherence time ~ 300 ms 
and hence a dynamic range of two oreders of magnitude. In the quadrupolar system described in 
this paper, the quadrupolar coupling value is ~ 16 kHz yielding a quadrupolar evolution gate time 
~ 20/is. Since the coherence times are 14ms and 4ms for inner and outer transitions respectively, 
the system yields a dynamic range of 3 orders of magnitude. Thus, both the dynamic range and the 
clock speed (inverse of gate time) are better in the quadrupolar system described here. 

V. CONCLUSIONS 

The implementation of quantum algorithms on quadrupolar nuclei validate their use as a alternate 
candidates for quantum information processing. DJ algorithm has been implemented here in a spin- 
3/2 system by manipulation of coherent superposition using evolution under quadrupolar interaction 
and r.f. pulses. The errors in our experiments were mainly caused by relaxation and imperfection of 
r.f. pulses. Use of tailored multi frequency pulses [22] can further decrease gate time and relaxation 
errors. Some quantum algorithms like Grover's search algorithm require uniform superposition of 
states, which can be realized by application of multiple quantum pulses [29]. Efforts are ongoing 



in our lab to develop such pulses and implement various quantum algorithms in 3/2 and 7/2 spin 
systems. Since completion of this work, an implementation of continuous version of Grover's search 
algorithm has also been reported in a spin-3/2 system [33]. ' 
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FIGURE CAPTIONS 



FIG. 1: (a) Energy level diagram of a spin-3/2 nucleus oriented in a liquid crystal matrix. The 
different spin states can be labeled as states of a 2-qubit system. The equilibrium deviation popu- 
lations of different states under high-field high-temperature approximation are schematically shown 
by the dots on the right hand side, (b) The equilibrium spectrum of 23 Na obtained after a hard 
( / k/2) v pulse. Along x-axis are the frequencies in KHz and along y-axis are the intensities. The three 
single quantum transitions are well separated by an effective quadrupolar coupling (A) of about 16 
kHz. The outer lines are broader than the inner line (line widths are different due to differences in 
relaxation matrix elements as well as due to fluctuations in S values). These fluctuations in S values 
affect only the outer transitions in the first order, reducing their transverse relaxation time T^. The 
integrated intensities are in the correct ratio of 3:4:3. The spectrum is plotted with a Lorentzian 
line-broadening factor of 200Hz. 

FIG. 2: Quantum circuit for implementing Deutsch-Jozsa algorithm. The first (ir/2)_ y hard pulse 
creates superposition of all states. Uf is the unitary transform corresponding to the function /. The 
last step is measurement. In NMR this step can be simplified to acquisition of signal immediately 
after Uf is implemented. The sign of the input qubit's resonances with respect to those of the work 
qubits resonances (Eq. [6]) distinguish between the constant and balanced functions. 

FIG. 3: Implementation of DJ algorithm on 23 Na (spin-3/2) nuclei. The algorithm starts from 
1 00) pseudopure state. After running through the quantum circuit of Fig.l, the acquired signal is 
Fourier transformed. The spectra corresponding to the operations U±, U2, U3 and C/4 are given. 
Along x-axis are the frequencies in KHz and along y-axis are the intensities. Constant functions are 
distinguished from balanced functions by the sign of resonance of input qubit. As seen by the single 
quantum coherences of Eq. [6]; for U\ and U 2 the resonance of input qubit (central transition) has 
same sign as the resonances of work qubit (outer transitions), implying that the corresponding func- 
tions f\ and f 2 are constant, whereas for U3 and C/4, the sign of the central transition are opposite, 
indicating that fs and f± are balanced. 
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